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In the present paper the model interaction potential between dust grains is proposed. It is pre- 
sented by two independent components, correspond to repulsive and attractive parts. The repulsive 
part represented by hard-core potential and the potential of Coulomb interaction. The attractive 
part is described by multipole interaction between two cells (we use the cell model of dusty plasma) 
which is long-range on its nature. The multipole interaction, by-turn, is presented by two addi- 
tive contributions: the interaction multipole — induced multipole, which is an analogue of dispersion 
interaction between neutral atoms (molecules), and the interaction between two multipoles. The 
dominating contribution to attractive part of interaction potential is determined. 
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I. INTRODUCTION 

Dusty (complex) plasma consist of weakly ionized gas (plasma) and highly charged microparticles (grains). There 
were several surprising effects discovered in complex plasmas: ordered structures and phase transitions in it, dust 
voids [1-3]. Therefore the problem of interaction between dust grains is significant. 

The forces, that order dust grains, are not correspond to direct Coulomb interaction (see, e.g., Ref [4]). Dust grains 
in a plasma due to various conditions gain huge charge (~ 10 4 in units of electron charge). This charge can be positive 
or negative. The main difficulty is to describe the attraction between between grains in such a system. 

There were a lot of papers, where this problem in one way or another was considered [5-9]. However, this question 
is still open. 

In the present work the model interparticle potential that predicts attraction between dust grains is proposed. The 
essential role in our construction plays the cell approach to dusty plasma. 

II. AVERAGED POTENTIAL OF MULTIPOLE INTERACTION BETWEEN DUST GRAINS 

In this section we formulate the general approach and discuss the method of its solution. 

A. The cell model of complex plasma 

According to its electroneutrality plasma can be represented as the collection of electroneutral cells of radius 

If 3 N 1/3 
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where n p is the dust grain density. The electroneutrality condition for a cell is 



(1) 



r c 

Ze + 4tt J p(r)r 2 dr = 0. (2) 

Here Z is theaverage dust grain charge, e is the electron charge, r p is dust grain radius, p(r) is the volume-charge 
density that can be determined in the self-consistent approximation. Note, that in the absence of an external electric 



'Electronic address: mishagliOonu.edu.ua 



2 



field the electric charge distribution has s spherical symmetry: p(r) => p(r). Eq. (2) allows us to establish dependence 
of Debye radius rpj on the average dust grain charge Z. This was discussed in detail in Ref. [10]. 

In the self-consistent field approximation the interaction between cells is absent. But it can appear due to fluctua- 
tions. Let us consider this mechanism in detail. 



B. Main statements of our approach 



Plasma represents as the collection of noninteracting cells in the self-consistent field approximation. The systematic 
perturbation of cell's spherical form occurs due to the thermal motion of dust grains and plasma particles (i.e. 
electrons and ions). Therefore fluctuation electric multipole moments occur and generate long-range electric fields. 

The electric field of one cell acts on its neighbours and generates interaction effects, which are similar to dispersion 
interaction between neutral atoms (molecules). In detail there are two interaction mechanisms: 

1. The fluctuation field of one cell polarizes neighbour one. Therefore, the second cell gains certain value of induced 
multipole moments M^- lnd \ Induced multipole moments M( ind ) interact with multipole moments M of the first 
cell. The average value of such an interaction is an analogue of the dispersion interaction between neutral atoms 
(molecules) : 

*1 - (WWnd>>. (3) 



2. The interaction of fluctuation multipole moments of both cells. The average value of such an interaction is 
determined in a different way: 

$2 = ~P(W 2 MM ,). (4) 
Here /3 = 1/kT, prime denotes multipole moments of the other cell. 

It is easy to ascertain, that both components lead to attraction between dust grains. The dominating contribution 
in both cases determines by dipole— dipole interactions, decreasing as 1/R 6 . At the same time, with the dust grains 
approach, when their "coats" of electrons begin to cover, repulsion forces are occur since grains have charges of equal 
signs. This allows us qualitatively to model the interaction between two cells (i.e. grains) with potential 

U(R) = U T (R) + U a (R), (5) 

where U r (R) and U a (R) are the repulsive and the attractive parts of U(R) respectively. The repulsive part of U(R) is 
modeled by the combination of the hard-core (the radius of hard-core coincide with grain's radius r p ) potential and 
the potential of Coulomb repulsion on distances 2r p < R < 2(r p + m)'- 

!oo, i?<2r p ; 
(Ze) 2 , , (6) 

The attractive part U a (R) is conditioned by forces of electric multipole interaction: 

U a (R) = $ 1 (R) + $ 2 (R). (7) 

The proposed interaction potential can lead to ordered structures for certain values of temperature and dust grain 
density. Let us determine dominating contributions to <f>i(i?) and <&2(-R) according to (3) and (4). 



C. "Dispersion" interparticle interaction 

The average energy of electrostatic multipole-induced multipole interaction is determined according to (3) as follows: 

$i (R) = (W DD(ind) ) + <W D Q ( .nd) ) + (W QQ(M) ) + ..., (8) 

where contributions (W DD (ind)}, (W D Q(ind)) and (VFggtind)) describes averaged dipole-dipole, dipole-quadrupole and 
quadrupole-quadrupole interactions respectively. 

The main contribution (W DD (ina) ) can be estimated from the following simple considerations. It was noted above 
that fluctuation field of one cell polarizes the other one. Therefore the second cell acquires induced dipole moment 
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d( ind ) = ad/i? 3 , where a is the polarizability of a cell and d is fluctuating dipole moment of the first cell. Thus the 
energy of dipole-induced dipole interaction is W DD ^&) = — d • E( ind ) (here E^ nd ' = d( ind )/i? 3 is the induced field of 
the second cell). The careful analysis of this problem gives 



~ -4a 



R e ' 



(9) 



D. Electrostatic multipole interaction 



The average value of the direct electrostatic multipole interaction given by (4) can be represented as follows 



(10) 



Here coefficients Ai (i = 6, 8, 10, . . .) express through the mean-square values of the multipole moments of a cell. To 
find dipole-dipole contribution we can follow the authors of Ref. [11] where dipole fluid is considered. Thus 



2\2 



<P 2 (R) ~ 



The rest contributions to (8) and (10) require more detailed analysis and will be considered in further paper. 



(11) 



III. POLARIZABILITY OF A CELL 



In the preceding section the model potential of interparticle interaction U(R) has been introduced. Also it was 
shown that we need to determine polarizability and dipole moment of a cell to know the main contribution to attractive 
part of U{R). Thus let us consider the electroneutral cell in an external electric field E . 

Under such conditions the electroneutral cell acquires polarization vector P = aEo. On the other hand, polarization 
vector is a dipole moment of the unit cell volume: 



= J r ■ p(r)dV, 



(12) 



where vector r is measured from the middle of the grain, p(r) is the volume charge density determined in the self- 
consistent field approximation: 



p(r, i?) = e 



n a o exp 



kT 



- rieo exp 



e</)(r, 1?) 
kT 



(13) 



Here z a e is the charge of an ion of kind a (e is the electron charge, z a are the positive and negative numbers), (f> is the 
electric potential inside the cell, n e o and n a o are the electron and ion of kind a densities respectively with = 0. It 
was taken into account that electric potential distribution in the presence of an external electric field loses its spherical 
symmetry. 

The potential (j>(r, #) satisfies the Poisson equation 



A</>(r,?9) = -47rp(r,t?), 
which is completed with the next boundary conditions: 



0(r c ,0) 



dn 



-Eqt c cos-d, 
— j + E cosi!). 



(14) 



(15) 



In further investigation we use the next dimensionless variables: 

r r c ~ eE 



o 



kT rp ' 



7 kT 



(16) 
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It is suggested that e<j> <C kT and the Poisson equation (14) can be linearized. Then the solution of (14) for 
renormalized dimensionless potential 



kT n c0 + ^ z^rieo 



(17) 



has the form 



Here 



V>(r,tf) = Vo(f) + Vi(M)- 



, , „ sinhf/A „ cosh fl A 



is the isotropic part of the potential, the angular part is proportional to cos $ and equals to 

tp(f, #) = E cost? 



^ 2 sinhf/A +A coshr7A 



, (| j A sinhf/A _ a2 cosh f /A 



(18) 



(19) 



(20) 



is the potential of an external field inside a cell and A = /r p is the dimensionless Debye length 



kT 



47re 2 ( n c o + E z l n aO 



rn = 

for the dust grain. Coefficients A i7 Bi (i=0, 1) are determined by the boundary conditions (15) and have the form 
(Z/Z ) cosh f + 7]^ (sinhi - A cosh i) (Z/Z ) sinh f + tj^ (-A sinh \ + cosh ±) 



A sinh + cosh 



A sinh + cosh s -^- 



and 



- 2? 3 sinh i - iAi+I ^ cosh l sinh ? + A cogh ? 

= So- 



il 1 — -En 



A(2<r - 2A 2 - 1) sinh s^l + (2A 2 (cr - 1) + ?) cosh ^1 

2ALhI ^ sinh i - 2<^ 3 cosh j - X sinh ^ + ? cosh f 
'A(2<r-2A 2 - l)sinh^i + (2A 2 (<r - 1) + 7) cosh s^i 



En A*, 



E B*. 



Thus founded electric potential (18) inside the polarized cell allows us, finally, to determine it's polarizability. 
Indeed, linearized Poisson equation gives us 

kT 1 

and according to (12) and P — aE we obtain 

a = -J {a\ [ - (<r 2 + 3A 2 ) sinh \ + 3A<r cosh i + (1 + 3A 2 ) sinh \ - 3A cosh i)] 
3 L A A A A 

+ Si [3A? sinh { - (c; 2 + 3A 2 ) cosh \ - 3A sinh \ + (1 + 3A 2 ) cosh \] }. (21) 
A A A A J 



The dependence of a/r 3 on r c /r p is shown in Fig. 1. 

Next, we turn to the definition of the multipole moments of the cell. 
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FIG. 1: Dimensionless polarizability a/r^ as the function of r c /r p . The parameters are: T = 2000 K, r p = 10 4 cm. Here 
considered ions of only one kind with z a = 1 and, therefore, ij = 1. 

IV. MEAN-SQUARE FLUCTUATION MULTIPOLE MOMENTS OF A CELL 

Let us consider the electroneutral cell with the perturbed surface shape. We assume that the fluctuations of the 
radius r c are small and can be expanded in a series of spherical harmonics: 

r c (0, if) = r<°> + 5r{0, <p), 5r(#, if) = r p £ lnm Y nm (ti, <p). (22) 

nm 

Here is the nonperturbed cell's radius, and dimensionless coefficients j nm that determine degree of deflection of 
the cell's surface shape. According to the thermodynamic theory of fluctuations (see Ref. [12]) 

kT 

(\ lnm \ 2 )^-—, (23) 
where W nm are the expansion coefficients for the increment of the electrostatic energy inside the cell: 

«W = ]TT4« m |7„ m | 2 . ( 24 ) 

nm 

To determine coefficients (|7„ m | 2 ) and multipole moments let us consider the evident expressions for the potential 
of fluctuation electric field. 

A. Fluctuation electric field 

Let and out be the electrostatic potentials inside and outside perturbed cell respectively. The potentials <f>i u 
and <j) ou t satisfy equations 

A0 in = -Airp, (25) 
A^out = (26) 
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and boundary conditions 



d(j)ir 



dn 



= — 47T(T, 



^out 



dn 

r— »oo 



(27) 



> 0. 



<9n 

0out(r) 

Note, that potentials 0i n and O ut depend on all spherical coordinates. The dependence on coordinate ip is absent 
above. 

The electric potential inside the cell is perturbed as well as the cell's surface shape. Therefore it can be 
represented in the dimensionless variables (16) and (17) as follows: 

ipin(f, (p) = ipo(f) + 5ip(r, ■d, ip). 

Here ipoi?) is the equilibrium electric potential inside the cell has been determined above (see Eq. (19)), and 
5tp(r, 1?, ip) = eS(j)(f, fl, <p)/kT is the potential caused by deviation of cell's surface from spherical one. Thus Eq. (25) 
after linearization (we also assume e<f)- m <C kT) disintegrates in two independent equations: 

A^ = and AStp = -j^H- 

The equation for ipo has already been solved in the preceding section. The fluctuation part Sip due to its smallness 
can be expanded in a series of spherical harmonics, and we obtain 



/ttA 



oo n 



<ty(M,¥>) = V2f E E W*,¥>) C mn I n+i (r/X)+D mn I_ n _ i (f/X) 



n—0 m— — n 



(28) 



where I n+ i (x) and I_ n _i(x) are the Modified Spherical Bessel functions of the first and second kind respectively. 
The Laplace equation (26) in the exterior has an obvious solution: 



0out(r,i?,y) = E ^TT 1 "" 1 ^'^)- 

nm 

It is known that the coefficients fj, nm are the required multipole moments. 



(29) 



B. Evident expression for multipole moments 

To determine coefficients C nm , D nm we will do the next. Expanding the first boundary condition in (27) in series 
near the ri^ and restricting ourself by contribution of the first order we obtain 



1 cty p (f) 



r p df 



5r(ti,<p) = 0. 



(30) 



The second condition in (27) under the 



dtp p (r) 



dr 



-Attct 



gives us 



—6i()(r,&,<p) 



0. 



(31) 



r=l 
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Equations (30), (31) allows us to determine coefficients C nm , D n 



-An 



Cnm — Tn 



sinh^/A cosh^/A 

-A 1 



Bo 



( sinh^/A ^cosh^/A 
V <r <r 2 



, l n+ l(l/X) + nXI nH (l/X) 
(^/ A )t , ^-777^ - J n+i(?A) 



-An 



'7_„_s(l/A)+nA/_ n _i(l/A) 
sinh<r/A cosh?/A 



/sinh^/A cosh^/A 
i>o A = 



I_ TI _ i (l/A) + nA/_ n _i(l/A) 
} 4 + 3(l/A) + nA4 + i(l/A) 



M/x) 



Now, we use the third condition from (27). We have: 

_ kT dxp p 



dn 

'0out 



er p dr 



+ 



¥T_ d&ip 
er p dr 



dn 



1 (n + l)fi mn (m) 

n+2 c n+2 '» ^' ^ ) ■ 



v ran 

From the equality of corresponding expressions before spherical harmonics we finally obtain 

kT 



Moo = ~r p — C 

e 



sinh^/A cosh^/A 
-A = 1 



B, 



kT <r 



n+2 



p e n + 1 



Qn'Yn 



where 



/„+|(c/A)+n-7„ + i(c/A) 



/ sinh^/A ^cosh<r/A 
(n > 0), 



} (?/A)+n-/_ n _i(?/A) 



We are interested in the case n > 0. So, the mean-square value of the multipole moments according to (33) is 

2 



(I Mr, 



2 ) = (V 1 



e n + 1 



Qn (|7« 



(32) 
(33) 

(34) 

i 

(35) 



Let us determine {\j n m\ 2 } according to (23). The electrostatic energy W c \ can be determined in a standard way: 
W * = -^ y [V0 in (r, 2 dr=i- / d( ^/ sinlMl? / {V[fo(r) + <ty(r,i?,¥>)]}Vdr. 



The dominating term of the last expression (when we put Sr and 5<j> equal to zero) is the energy of the 

electrostatic field of charged dust grain obtained and discussed in Ref. []. Now, we are interested in the term (24) 
which corresponds to the fluctuation term 8<j>. To determine SW we take into account the smallness of Sr. Let us 
expand the last integral over the upper limit. Note, that contributions proportionate to Sr can be omitted due to 
(Sr) = 0. Thus, constraining with the first infinitesimal order and taking into account (24) we obtain: 

W nm = ^ (^) 2 r P (V« + W<» + W<&) , (36) 



where 



W (1) 

vr nm 

w {2) 

* v nm 



dip d 2 ^ ~ 2 

r 

dr dr 2 



dn I r 2 df I V 



dr 



dnY nm (0,<p)Y lp {#,<p), 



■^Y nm (0,<p) C n I n+h {f/\) + D n I_ n Af/X) 



= 2 



dipp d 
dr df 



C n I n+h {f/X) + D n I_ n _ h {f/X) 



f 2 
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Here V is an usual nabla operator V written in dimensionless form. 

C. The mean square value of dipole moment of a cell 

Let us put n = 1 and m = 0. Then 



and we have for (i — 1, 2, 3): 



(Z/Z ) 2 , z/z 



+ 



c£> 2 <^X 2 D 



Yio = a/— cost?, 

47T 



A(2^ - ^ 2 - 2A 2 ) sinh i-1 + (2U 2 - 2A 2 - c 2 ) cosh i-i 



TV 



+ 2^{ ^ + 1} + + A2) ^ - ^ Sinh2 V 

+ A [2 ? 3 + 4a 2 - (2? + A 2 )(A + 1)] cosh2^ + A 3 (A 2 - 1) J, 

= ^ Qc 3 - \? e\ 2 + ^ 3 A 2 + A 4 - , 3 A 4 ) (C 2 - £ 2 ) 

-<;Afc 2 + 8A 2 ) sinh2<rA - ^ 3 A(8A 2 + 1) sinh % + A 2 (^ 2 + A 2 ) cosh 
4 4 A A 



+ <r 3 A 2 (A 2 + l)cosh 



A 



(Cl + D\) 



TV, 



(3) 



1 



10 



(Z/Z )X - (? - 1)A cosh i—i + (A 2 - sinh i— ^ 

A A 



(<T 2 + 2A 2 ) sinh ^ - 2<rA cosh ^ (Ci + £>i) 

A 



A 

D = X sinh — - h cosh — - — . 



V. CONCLUSION 

In the work a new approach to description of interaction between dust grains in plasma is proposed. It is taken 
into account that grains are in the electroneutral cells in equilibrium. Therefore the direct interaction between grains 
is absent. However, a thermal motion of dust grains plasma particles leads to charge fluctuation inside a cell and to 
surface shape fluctuation of a cell. The latter type of fluctuations is discussed in the presented work. 

It is shown, that the deviation of the surface of a cell from the spherical form is accompanied with the occurrence 
of non-zero multipole electrostatic moments that form long-range fields. The interaction of fluctuation multipole 
moments of the neighbouring cells with these fields leads to interaction between the cells (i.d. the grains), the 
dominating contribution of which is decreasing as 1/i? 6 . At the same time it is clear that at the short distances 
between the dust grains the screening effects sufficiently decay and the interaction has the character of Coulomb 
repulsion of the two charged grains. 

The evident form of the potential of interaction between the dust grains is proposed. In principle such a potential 
can describe the ordering effects of the dust grains, in particular the formation of quasi-crystal structures. The detail 
analysis of the proposed potential will be the subject of further work. 
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